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Abstract 

In this communication we deal with the exact solutions called "pseudosimilarity" of a steady free convection prob- 
lem studied by by Kumaran and Pop (2006). They showed that there is no similarity solution for the case of a wall 
temperature as T^{x) ^ x~i (resp. a wall heat flux as qw{x) ^ x~i , and a dimensionless heat transfer coefficient 
hw{x) ~ a;^^). We shall present some results about existence and asymptotic behaviour of new exact solutions of the 
resulting boundary value problem for each case. 
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1. Introduction 

Heat and mass transfer in saturated porous media occurs widely in natural phenomena and indus- 
trial applications, such as extrusion of polymers, continuous casting, colling of metallic plates,etc. 
More details about the existing literature on this subject can be found in the books (Neild and Be- 
jan, 1999; Vafai, 2000; Pop and Ingham, 2001) and the references therein. 

Recently Kumaran and Pop (2006) presented some original results about a particular case of this 
phenomena. They considered a fluid (water) characterized by a relationship between the tempera- 
ture T and the density p 

^^ = -i{T-T,r, (1.1) 

pc 

where is the maximum density in the liquid phase, and 7 is the fluid thermal expansion coeffi- 
cient of water at 4°C. We note that equation (11.11) have been obtained by Goren (1966). This kind 
of problems have been the subject of many papers, see for instance (Soundalgekar, 1980; Black el 
al., 1984; Poulikakos, 1984). 

In (Kumaran and Pop, 2006), the authors studied the free convection about an impermeable flat 
plate embedded in a porous medium filled with water near the vicinity of its density maximum 
associated with the temperature of 4°C at atmospheric pressure with the condition u{x, 00) = 
(uniform free stream velocity). They showed analytically and numerically, the existence of multi- 
ple similarity solutions to the governing boundary layer equations for a class of problems namely, 
wall temperature (VWT, T^(a;) = x™), variable heat flux (VHF, Qwix) = x"^^), or variable heat 
transfer coefficient (VHTC, /iu;(x) = x^^). They also gave a nonexistence result for (m = — i), 
this singular case will be the main goal of the present investigation. We shall adopt a pseudosimi- 
larity approach to construct exact solutions for the resulting boundary value problem and to study 
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their precise asymptotic behaviour. 

According to Pop and Ingham (2001), the physical model is described by the following equations 

Ux + Vy = 0, 

u = T\ 

< UTr,+VTy =Tyy, (1.2) 

v{x, 0) = VwX^, 

u{x, oo) = T(x, oo) = 0. 
Supplemented by one of the following conditions 



T(x,0) = x™ {VWT), 

Ty{x,0) = x^ {VHF), (1.3) 

-^Ty{x, 0) = (VHTC). 



Where u and v are the velocity components along x and y axes, respectively, T is the fluid temper- 
ature and the exponent m is a real number. 

In the usual way, we introduce the stream function ip, (satisfying u = ipy and v = —il>x)- Then the 
above model can be expressed in a simpler form 

i^yTx - IPxTy = Tyy, (1.4) 

i^xix, 0) = -v^xP, limy^oo i^y{x,y) = 0. 
Define the similarity transformations 

i){x,y) = x"" firi), T{x,y) = x'^eir]), r] = yx\ (1.5) 

where / is the dimensionless stream-function, 6 is the dimensionless temperature, rj is the similar- 
ity variable, and the real numbers r and s satisfy the relations r + s = l,r — s = 2m and p = ^'"^""^ . 
Injecting (11.51) in (|1.4I) we get 

f = e\ 

9" + ^fe' -mfe = 0, (1.6) 
/(O) = 9{oo) = 0, 

subject to one of the following conditions 

^^(0) = 1 (a), 

^'(0) = -1 (6), (1.7) 

9'{0) + 9{0) = (c). 

Where = ~ 2m+i '^«> (th^ suction/injection parameter), and the prime denotes differentiation 
with respect to 7]. 

The boundary value problem (11.61) , for m = — |, under the condition (i), i = a,b, c, will be the 
studied in detail in the next section. 



2. The pseudosimilarity solution 

Magyari et al. (2002) have studied the well-known free convection boundary layer flow problem 
about a vertical flat plate with an inverse-linear temperature distribution T^(x) = Too + f where 
a > 0. under the condition u(x. oo) = 0. Thev have shown that the similaritv solutions in the 



classical fomi are missing, and proved that in order to overcome this difficulty, the stream function 
if) has to be shifted by the term a log(x), cr > in order to get exact solutions. In this section we 
adopt the same concept to obtain exact solutions for problem (|1.5I) - (|1.6I) . under the condition (11.41) . 
when m = — i. First, we study the case of a Variable Wall Temperature. 



• The (VWT) flow 



In this situation, the wall temperature is given by 

T^(a;)=x-i (2.1) 
Following the work (Magyari et al., 2002), we define the new stream-function 

ilj^x.y) = F{ri,x). (2.2) 

If m = — i we deduce from the previous section that(r = 0, s = 1) and then the pseudosimilarity 
transformations are defined by 

-^{x,y) = F{ri,x), T{x,y) = x~^<d{r],x), t] = yx~^ , (2.3) 

where F is the modified dimensionless stream-function, 6 is the dimensionless temperature and r/ 
is the pseudo-similarity variable, substituting in (11.51 ) we obtain 

F, = e\ 

(2.4) 

®W ~l~ 2^*?® -^Fri^x •^FxQ-q- 

As in (Ibrahim and Hassanien, 2001), the term 6^^^ can be neglected and we may assume that 
9 = Q{ri,x) (then the term F^^i can also be neglected). By writing F{r],x) = f{r]) + h{x), 
equation (11.41) yields 

6" + \f'e = -xh'ff. ^ ^ 

Hence there exists a real number a such that 

/' = e\ 

xh' = a, (2.6) 
0" + \f'e + aO' = 0. 

The integration of the second equation gives h(x) = a log(x) + est, and then the new stream- 
function is given by ^(x, y) = f{i]) + a log(x) + est. Coming back to problem (11.41) . we see that 
a has a physical meaning since a = —v^ (the suction parameter). 
In the remainder, we deal with the following boundary value problem 

/' = e', 

0" + \f'0 + ad' = 0, (2.7) 
9{0) = 1, e{oo)=0. 

Let us notice that the no-slip condition /(O) is not required, in fact it can be any real number (see 
(Hammouch, 2006) for more details). On the other hand, the substitution of (|2.7l) i in (12.71) 9 gives 
an unforced Duffing equation (Panayotoukanos et al., 2002) 

e" + -e^ + ae' = (2.8) 



with the conditions 



^(0) = 1, ^(oo) 



0. 



(2.9) 



whose study in the phase plane reveals the stability of the unique equilibrium point (0, 0) is related 
to the sign of a. Fisrt of all we show that a has to be positive, for this sake we define the Lyapunov 
function for equation (12.81) by 

. 1 



4 



-0^ 



(2.10) 



which is positive and satisfies V'ij]) = —aO'"^. If we assume that a < 0,V is monotonic increasing 
and then condition 9{oo) = could not be satisfied. Then, the real number a is assumed to be 
positive, it follows that the equilibrium point (0, 0) is globally asymptotically stable and all local 
solutions of (|2.8I) - (|2.9I) remain bounded and tend to zero at infinity. 

However, we remark that we cannot exclude the existence of oscillating solutions for smallest 
values of a (see Figure. 1). 




20 
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Figure. 1 Temperature Profiles as a function of rj for small values of cr and 0'{0) = —0.1, for the case (VWT). 



To establish the existence of monotonic decreasing positive solutions to (I2.8I) - (I2.9I) . we shall adopt 
a Crocco variables approach (Nachman and Callegari, 1980; Ishimura and Ushijima, 2004). For 
the sake of simplicity, we assume that the function 9{r]) (which can be regarded as an independent 
variable) is strictly monotonic (9' < 0, for 77 > 0). Actually, solutions which are increasing on 
some interval (0, r/o) may exist, although this case can be studied in an analogous way. 
Let r](9) denote the inverse function of 9{r]). Defining 



9, 



and 



leads to 



under the conditions 



Integrating over (0, s) leads to 



m = 



- -s^ < s < 1, 
and 0(00) = 0. 



1 r r3 , , , 

[s) = as-- / -T-:dr, s G (0,1). 



(2.11) 

(2.12) 
(2.13) 

(2.14) 



To classify the pseudosimilarity solutions to (|1.5I) - (|1.6I) . we study the solutions of the more simpler 
integral equation (12.141 ) with the conditions (12.131) . We notice that there is some analogy between 
problem (|2.8I )- (|2.9I ) and the following one: 

Find the traveling-wave (TW) solutions U (x, t) = 9(x — at), with the condition 9(oo) = 0, (where 
a is the wave speed), for a special case of the Fisher-KPP equation 

Ut = U^, + U\ (2.15) 

In the light of the work by Gilding and Kersner (2005) concerning this kind of problems (2.15), 

we deduce the existence of a threshold value a* such that existence and uniqueness of solution to 

(I2.8I) - (|2.9I) is guaranteed if cr = a*. While for a > a* there are multiple solutions. 

Now we pay attention to get estimation for a*, and to study the precise asymptotic behavior of 9 

, . 2 

as — > oo. Multiplying equation (12.121) by - leads to 

s 

02 y (0_^s)2 



then 



— 1 +- ^ + ^2 = ^2, 

si 



— + - < a'. 



A simple integration over (0, 1) gives 



consequently 



4,{1)' + ^<<t\ (2.16) 

Now we proceed to determine the large ?7-behavior of these solutions. For this sake, we will exploit 
an idea of Brezis et al. (1986). First we show that for large r] we have the following limits 

9' 9' 
lim 77 = -a (^) or lim - = (^). (2.18) 

rj—too 9 Tj—too 9 



Actually, setting xiv) = gives 



x' + ^X + x' + y = 0, (2.19) 



for all r] E (0, oo). Because that 0(s) < as we get that the function x is negative and bounded 
(.—o' < xiv) < 0)' thsii for V large x has a finite limit, say /. This limit satisfies P + la = which 
leads to (12181) . 

Assume now that (*) holds, then we have for all s G (0, 1) 

a] = - -^Q ^ ' ^ " . (2.20) 



s 

By the L'Hopital rule, we get 

(2.21) 

s^o s'^ s ba 



then 

(p(s) ~ as + — as s ^ (2.22) 



s' 



Next, we assume that (i^ic) holds, from equation (12.121) we get immediately 



as s ^ 0. (2.23) 

2cr 



Coming back to problem (|2.8I) - (|2.9I) we conclude that 

- For a = a* : 6 e ^ (exponential decay). 

- For a>a*: (algebraic decay). 



In the remainder we shall obtain some estimates for 0(1). We stress that the real number = 
—9'{0) plays the role of the shooting parameter for problem (I2.8I) - (|2.9I ). We note also that equation 

(12.161) gives a lower estimation for 6'{0), actually we have 6'{0) = —0(1) > — -y^cr^ — -. 

To obtain an upper estimation for 9'(0) we look for a positive solution to (|2.8I )- (|2.9I ) such that 
( := 9' (0) > 0. In the above analysis, we have said that for any C 7^ the (unique) local solution 
to (I2.8I) - (I2.9I) with the initial condition 6(0) = 1, is global and goes to zero at infinity. Since ( 
is positive there exists a positive number 77^ such that 9 is monotonic increasing on (0, 77^) and 
^'(^c) ~ 0. Setting A = 9{r](^) and 9(^ = 9(r] + 77^) for all i] > 0. The shifted function 6*^ is a 
solution to (|2.8I) - (|2.9I) with the initial condition 

e^{o) = A > 1, ^'^(0) = 0. 

Using the transformations 

T] — > A^r/, and ^ ~f ' 

A 

Clearly the function 6^^ is a solution to the following problem 
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9'i + cTx9'^ + y-^ = o, 

(2.24) 

^^(0) = 1, ^^(0) = 0, 9^{oo) = 0, 
where ax = A^^cr. 

Thanks to the above analysis, problem (|2.24|) has positive solution if and only if 

. > ^. ("s 

On the other hand, integrating (12.241) over (0, r^^) gives 

From which we deduce that aX < C + a. Hence inequality (|2.25l) is satisfied if 

C < a (aVs - 1) . (2.26) 
In conclusion, under the (VWT) condition, problem (|2.8l )- (|2.9l )-(a) admits multiple solutions (flows). 



Every solution is uniquely parametrized by the skin-friction parameter ^' (0) G — y cr^ — | , cr ((t-\/3 — l) 
The above results have been illustrated numerically in Figure. 2. 



1.2 
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Figure. 2 Temperature Profiles as a function of rj for various values of a and related values of 9'{0), for the case 

(VWT). 

• The (VHF) flow 

This flow is characterized by the wall heat flow = and the wall temperature T^(a;) = 0. 

In this case, we study problem (|2.8I) . (|2.9I) supplemented by condition (b). 

First we note that if 6^ is a solution to (I2.8I) - (|2.9I) then (—9) is also a solution. Hence The above 
analysis for (VWT) can be extended to the case of a (VHF) flow. Consequently, multiple solutions 
exist if a lateral injection is applied with an injection parameter sufficiently large (v^ = a > ■^). 

• The (VHTC) flow 

Consider now the case of (VHTC) flow (|2.8l) - (|2.9l) -(c). In such situation the dimensionless heat 
transfer coefficient is given by h^^^x) = x^^. Taking into account condition (c), we look for 
solution satisfying 

^(0) = -^(0) = K, 

where n is positive constant. Introducing the new variables 

0(r}) 

rj K^T], and 6 

n 

we find the new problem 

^" + £ + ^' = with = 

(2.27) 

^(0) = 1, ^'(0) = e{oo) = 

Using the same arguments as for the (VWT) flow, we deduce that there exists a minimal value a* 
such that problem (12.271) has positive solutions only if > cr^*, more precisely a* > -j^. Every 



solution is parametrized by ^^'(0) G — y o"k ~ (o"kV^ — l) 



3. Conclusion 



In the present paper, we have showed how to construct exact solutions by introducing an additional 
logarithmic term in the usual stream-function. The relevant problem has been studied via a Crocco 
transformation combined with an integral equation method. The following conclusions have been 
made as a result of our investigation: 

• There is existence of a minimal value a* (suction or injection) such that multiple solutions exit 
only for a > a* and uniqueness holds for a = a*. 

• Solutions for a — a* exhibit an exponential decay while for a > a* all solution are decaying 
algebraically. 
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